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ABSTRACT 

We present a method for computing the net transmission of X-rays emitted by shock-heated plasma 
distributed throughout a partiahy opticahy thick stehar wind from a massive star. We find the 
transmission by an exact integration of the formal solution, assuming that the emitting plasma and 
absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which 
there is assumed to be no emission. This model is more realistic than either the slab absorption 
associated with a corona at the base of the wind or the exospheric approximation that assumes that 
all observed X-rays are emitted without attenuation from above the radius of optical depth unity. 
Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined 
table of the wavelength dependent wind opacity. We provide a default wind opacity model that is 
more representative of real wind opacities than the commonly used neutral interstellar medium (ISM) 
tabulation. Preliminary modeling of Chandra grating data indicates that the X-ray hardness trend of 
OB stars with spectral subtype can largely be understood as a wind absorption effect. 
Subject headings: radiative transfer — stars: early type — stars: mass-loss — stars: winds, outfiows 



1. INTRODUCTION 

The absorption of soft X-rays by the powerful, 
radiation-driven winds of OB stars has long been rec- 
ognized as a significant effect both on the X-rays ob- 
served from these stars and on the physical condi- 
tions in their winds. The soft X-ray emission ob- 
served in OB stars by Einstein and ROSAT implied 
only modest wind attenuation of the X-rays, and thus 
ruled out significant coronal emission as a source of 
the ubiquitous X-ray emission seen in these massive 
stars (Cassinelh & Oison 1979; Cassinelh fc Swank|[T98l 
Eillier et al. 1993; MacFarlane et al. 1994). For this and 
other reasons, the wind-shock paradigm for the produc- 
tion of X-rays in O B stars has become general ly ac- 
cepted (e.g.. .Owocki et al.) 119881: iPauldrachet Id] 119941 : 
iFeldmeier et al.l Il997allbl : IXahn et al.l I2QQ1D . "although 
many aspects are still poorly understood. 

Using X-ray observations to constrain X-ray produc- 
tion mechanisms requires proper account of the ab- 
sorption of distributed sources of X-ray emission. The 
amount and wavelength dependence of the wind ab- 
sorption can be used as a diagnostic of the loca- 
tion/distribution of the shock-heated plasma and of the 
wind mass-loss rate, especially in terms of its effect on in- 
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dividual line profile shapes (|Owocki fc Cohenll2QQlL here- 
after OCOl). Even simply deriving an intrinsic X-ray lu- 
minosity for energy budget considerations requires cor- 
rectly modeling the significant attenuation of the emit- 
ting X -rays, especially in the dense winds of O super- 
giants (|Hillier et al.lll993l : lOwocki fc Cohenlll999[ ). 

Because the emitting plasma is spatially distributed 
throughout the absorbing wind, simple prescriptions for 
the attenuation can be inaccurate. Specifically, the com- 
monly used slab model of absorption, appropriate for an 
intervening interstellar medium (ISM) cloud where all of 
the emission originates beyond the absorbing medium, 
has transmission with an exponential dependence on slab 
optical depth, and strongly overestimates the amount of 
attenuation as the wind becomes optically thick. Fur- 
thermore, because hydrogen and helium are ionized, 
models for neutral gas significantly overestimate wind 
opacity. Unfortunately, due to the lack of easily available 
tools for incorporating appropriate radiative transfer and 
opacities, inadequate models, such as those intended for 
neutral interstellar medium absorption, are routinely ap- 
plied to account for wind attenuation. This is in spite of 
the fact that a number of previous works have recognized 
the necessity for an d applied appropr i ate wind absorp- 
tion prescriptions (jHillier et aTl 119931: iPauldrach et all 



l2QMI : 1nwQcki fc Cohenll2QQll : lOskinova et al.l[2QQ6h . 

We have developed a method for implementing an ex- 
act solution to a realistic model of the radiation transport 
that can be easily combined with a pre-calculated opacity 
table to find the wavelength-dependent emergent X-ray 
fiux from a stellar wind. We also provide a reasonable 
default opacity model that can be used for most OB star 
winds. Our analysis tool, which we name windtabs^ for 
wind ^able a^^orption, can be combined with an indepen- 
dent plasma emission model, su ch as the As trophysical 
Plasma Emission Code (apec) (jSmith et al.l 1200 1 ) that 
is widely used in fitting stellar X-ray spectra. This can be 
used to realistically model the low-resolution CCD spec- 
tra that are produced in large quantities by surveys of 
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clusters and OB associations with Chandra and XMM- 
Newton fe.g. JSana et alJl2QQ6l : IWang et alJl2QQ8D . It can 
also be used to model grating spectra in detail, and pro- 
vides a means of disentangling the wind absorption ef- 
fects from the emission temperature effects that both ap- 
pear to contribute to the recently discovered trend in the 
morphology of OB star spectra observed at hi^h reso lu- 
tion with the Chandra gratings (|Walborn et al.ll2QQ9l ). 

2. RADIATION TRANSPORT MODEL 

In this section, we derive an expression for the frac- 
tion of X-rays transmitted from a massive star wind 
as a function of the wavelength-dependent opacity, the 
mass- loss rate, and the wind velocit y law. We make 
assumptions similar to those made in lOwocki fc Coheiil 
([2001): we model the wind as a spherically symmetric 
two-component medium, where a small fraction of the 
wind is heated to X-ray emitting temperatures (Tx ~ 
1 — 10 MK), while the bulk of the wind is composed of 
relatively cool material (Twind/^eff ^ 0-5 — 1) that can 
absorb the X-rays via the bound- free opacity of the mod- 
erately ionized metals. We assume that the X-ray emis- 
sion turns on at some radius Rq > R^. We also assume 
that the temperature distribution of the X-ray emitting 
plasma is the same over the entire emitting volume. 

The assumptions we make regarding the distribution 
of the X-ray emitting plasma are rooted in the available 
observational evidence, as well as in the results of ex- 
tensive theoretical simulations. Detailed studies of emis- 
sion line profiles as well as constraints from forbidden- 
to-intercombination line ratios in He-like ions support a 
picture in which X-ray emission starts at Ro ^ l.bR^^ for 
all observable ions ( Cohen et al. 20061: [Leutenegger et al.l 
|2006"; 'Co hen et al.l 1201 0). with a roughly constant fill- 
ing factor above the onset radius. Numerical simulations 
have typically found onset radii for stro n g shocks that are 
compar able to this (.Owocki et al.lllQSSl ). iFeldmeier et al.l 
(|1997b[ ) performed simulations of winds seeded with 
base perturbations. They have found that clump-clump 
collisions are important, and that the resulting X-ray 
emission is distributed pri marily below 10 stellar radii. 
iRunacres fc Owockil (|2002l ) have also found X-ray emis- 
sion with onset radii of order 1.5 stellar radii and extend- 
ing out to large radii. However, the results for the onset 
radius are also sensiti ve to the treatment of the scat tered 
radiation field (O wocki fc Pulslll996l : [Owoc ki"2009l. 

The observed X-ray luminosity as a function of wave- 
length is given by 
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where r]x{r) is the X-ray emissivity, and r{r^fi,X) is the 
continuum optical depth of the dominant cool component 
along a ray from the emitting volume element to the 
observer. 

The optical depth can be derived as in 
lOwocki fc CohenI (|2001[ ) for a smooth wind. It is 
given by the integral 



r{p,z,X) 
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Here p and z are ray coordinates, with impact parameter 
p = Y^l — jj? r and distance along the ray z = jir^ where 



ji is the direction cosine to the observer at local radius 
r. n{\) is the atomic opacity of the wind, and p{r) is 
the density of the wind. Using the continuity equation, 
p{r) = M/4:7Tr'^v{r), where M is the stellar mass loss 
rate, and defining the characteristic wind optical depths 



we can write 
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is a dimensionless integral that depends purely on the 
ray geometry. Here, R^ is the stellar radius, ^oo is the 
wind terminal velocity, and w{r) = v{r)/voo is the scaled 
wind velocity. Note that we have assumed that n{\) is 
constant throughout the wind; we further discuss this 
assumption in Section [3l We take the velocity to follow 
a beta law: v = i'oo(l — R^/r)^. We also take /3 = 1 
in this paper as a good approximation for many O star 
winds; however, evaluation for general values of P is not 
difficult. 

The emissiv ity is assumed t o sca le with density 
squared, as in lOwocki fc CohenI ()2001f ). We ignore the 
Doppler shift of the emitted X-rays and write 



r]x{r > Ro) = r]x,o 
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This expression gives the correct radial dependence for 
the total line strength, but discards information about 
the profile shape. This is a justified approximation in cal- 
culating the broadband X-ray transmission of the wind. 
Here we assume that X-ray emission begins at a mini- 
mum radius Rq^ with po = p{Ro) and 7?a,o = Vxi^o)- In 
this paper, we will assume that the X-ray filling factor 
is constant with radius; it is trivial to add a power-law 
radial dependence, as in OCOl. We also assume that X- 
ray emissivity follows the same radial distribution at all 
observable wavelengths. 

The model described in the preceding paragraphs is 
illustrated in Figures [1] and [2l The model emissivity and 
transmission are visualized separately in Figure [U and 
together in the left panel of Figure [2j The right panel 
of Figure [2] gives the net transmission for the exospheric 
approximation for comparison, which we discuss at more 
length below. 

Using Equations ([I])-([6]), we can calculate the trans- 
mission of the wind as a function of r* . The transmission 
is simply the observed flux (Equation ([1])) divided by the 
unattenuated flux, which can be found by setting r* = 
in the same equation: 
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To account for the attenuation in the numerator, it is 
convenient to define an angle- averaged transmission from 
each radius r. 



1 f^ 
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is thus 
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Figure 1. This diagram shows the X-ray emission and absorption 
properties of a model stellar wind. The observer is on the left. The 
black disk denotes the star, and the thin black line shows the onset 
radius for X-ray emission (1.5 i?* in this case). Contours of con- 
stant X-ray emissivity (proportional to density squared) are shown 
with dashed lines at intervals of an order of magnitude in differen- 
tial emissivity. The absolute scale is normalized to the maximum 
emissivity at the onset radius of X-ray emission. Contours of con- 
stant continuum optical depth calculated for r* = 3 are shown with 
solid lines. 
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Figure 2. Left panel is the same as Figure \T\ but plotting the 
net emissivity (product of intrinsic emissivity times transmission, 
including occultation by the stellar core) instead of the intrinsic 
emissivity; this net emissivity has been renormalized so that the 
point of maximum net emissivity has a value of unity. The right 
panel shows the same plot in the exospheric approximation, where 
the gray zone is inside the radius of (radial) optical depth unity 
and thus does not contribute to the observed X-ray flux. 



where 



M* 



1 



m 



(9) 



gives the ja coordinate of occultation by the the stellar 
core. Some X-rays are obscured by the stellar core even 
when the wind is transparent: 



To(r)=T(r,0) = -(l+M*)- 
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We can further evaluate this expression by substituting 
the continuity equation, and by defining the inverse ra- 
dial coordinate u = R^/r: 
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Figure [3] shows the angle- averaged transmission T{u) 
for selected values of r* . Figure S] shows the following 
quantity: 



Lx{u,T,) _ j;du'w-^{u')T{u',n) 
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This can be thought of as the cumulative distribution 
of observed X-ray emission, starting from u = {R = 
oo) and integrating in to u = uq. Together, Figures [3] 
and [4] show the relative importance of transmission and 
emission as a function of radius in stellar winds. For the 
entire range of interest in r* , essentially all of the wind 
down to uq contributes to the observed X-ray flux. 

Figure [5] shows T{r^) for this model, along with com- 
parisons to two other absorption prescriptions: a simple 
intervening absorber, T = e~^ , appropriate for a coronal 
slab model, e.g., as implemented in the XSPEC mod- 
els wabs or tbabs 'dind a n exospheric approximation (e.g. 
lOwocki fc Cohenlll999l ). where T = below the radius 
of optical depth unity, and T = 1 everywhere above it: 
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where Ux = m.m{ui{r^)^uo). The inverse radial coordi- 
nate of optical depth unity is given by evaluating the 
optical depth integral (Equation (J5])) along a radial ray 
(p = 0, z = r), with the result 



^i(r*) = 1 



-l/r* 



(15) 



for /3 = 1 . Note that in Figures [SJISI we have used /3 = 1 
and Ro = 1.5i?*. We stress that Equation (fT4|) is simply 
the consequence of using a step function for the angle- 
averaged transmission T in Equation (p!2)) . rather than 
the more realistic expression given in Equation (|8]). 

The transmission of our model falls off much more 
gradually than e~^ (Figure [5j), but it is also more ac- 
curate than the exospheric approximation, especially at 
moderate optical depth. The exospheric approximation 
has the correct asymptotic behavior for large values of r* , 
but overestimates the transmission by a flxed factor. The 
exospheric transmission at large r* can be brought into 
agreement with windtabs by multiplying the exospheric 
r* by three. 

It would be possible to generalize the radiation trans- 
port model described in this section to include poros- 
ity by introd ucing an appropriate defl nition of effective 
opacit y, as in lOskinova et al.l (|2QQ6f ) or lOwocki fc CohenI 
However, a speciflc implementation of this and 



Integrating over shells at all radii, the net transmission discussion of its consequences is left to a future study. 
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Figure 3. Angle averaged transmission T as a function of shell 
inverse radial coordinate u = R^/r^ with different curves for dif- 
ferent characteristic wind optical depths r*. The equivalent plot 
in the exospheric approximation would be a step function at wi, 
the inverse radial coordinate of optical depth unity. Note that the 
transmission is less than unity even for r* = due to occultation 
by the stellar core. 
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Figure 4. Fraction of X-ray emission originating from outside 
inverse radius w, normalized to total X-ray emission integrated to 
no = 2/3. For all but the very optically thick case of r* = 10, the 
emission comes from a wide range in u. 

3. OPACITY MODEL 

To determine the wind transmission as a function of 
wavelength, we need to account for the wavelength de- 
pendence of the optical depth through the opacity, which 
we write as 

r*(A) =/>:(A)I]*, (16) 



where 



M 



AttR^Vo 



(17) 



is the characteristic mass column density of the wind (in 
gcm~^). 

The continuum opacity of a stellar wind in the X-ray 
band can be calculated by summing the contributions of 
each constituent species. Thus, we must know the atomic 



Figure 5. Comparison of transmission of three different mod- 
els: coronal slab (e~'^), exospheric, and more realistic wind model 
(windtabs) . The fixed parameters are ^ = 1, i?o = 1-5. 

cross-sections, ionization fractions, and elemental abun- 
dances. Of these three, the cross sections are known 
to sufficient accuracy (e.g. Vern er fc Yakovlevlll995f ): the 
ionization balance may contribute some uncertainty in 
the calculation of the opacity, but is usually not the dom- 
inant source of error; and uncertainties in the elemental 
abundances are typically the most important. 

The opacity due to photoionization of a given shell of 
any individual species scales approximately as Hii{X) oc A^ 
above the threshold energy ^th of the shell. Because mul- 
tiple species are usually important for the X-ray opac- 
ity of astrophysical gas, the run of opacity with wave- 
length has a characteristic sawtooth shape, with individ- 
ual teeth at the ionization threshold energies of dominant 
ionization stages of abundant elements. 

O star winds are photoionized, with H and He fully 
stripped (although He may recombine in some dense 
winds; see below), and most other elements mainly in 
charge states +3 and +4. Thus, the opacity of stellar 
winds in the range 1 A < A < 40 A is dominated by 
K-shell absorption in C, N, and O, since they are the 
most abundant elements. Significant contributions from 
K-shell absorption in Ne, Mg, and Si are also present, as 
well as Fe L-shell absorption. 

The opacity of adjacent ionization stages of the same 
element are usually comparable, with the exception that 
the photoionization threshold energy is shifted. The ef- 
fect of a moderate shift in ionization balance on the opac- 
ity is relatively minor, since none of the strong X-ray 
emission lines observed from O star winds falls between 
the threshold energies of adjacent stages of abundant 
ions. However, the difference in threshold energies, and 
hence the broadband absorption, between neutral and 
wind material is significant. 

Therefore, while it is important to use an appropriate 
model for the wind ionization, it is sufficient for many ap- 
plications to use a single approximate ionization balance 
to model all O star winds. This is true even though the 
ionization balance can vary to some extent with radius, 
and also is different in different stars. 

If it is desired to model the opacity of a particular 
star, it is possible to construct a detailed model opac- 
ity for a stellar wind by using the output of a radia- 
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tive transfer model such as CMFGEN (|HimeLela]J [1993|; 
iHillier fc MiHerlfToOSl : iCohen et alJ[2QlQL J. Zsargo et al., 
in preparation). 

To illustrate the importance of various assumptions in 
opacity modeling, Figures [6] and [7] compare several model 
opacities. Figure [6] compares neutral interstellar medium 
opacity and a simple model for the op acity of a stel- 
lar w ind. Both assume solar abundances (jAsplund et al.l 
|2QQ9[ ). In the wind model we assume an ionization bal- 
ance with hydrogen and helium fully stripped, oxygen 
and nitrogen in the +3 charge state, and all other ele- 
ments in the +4 charge state. This is a good approxima- 
tion to the ionization balance in a typical O star wind, 
which is set by photoionization from the photospheric 
UV field. The model wind opacity is much lower than 
the model ISM opacity, especially at long wavelengths, 
mainly due to the ionization of hydrogen and helium. 
The shift in ionization threshold energies due to the pres- 
ence of more highly ionized species is also clear. 

Figure [7] shows three stellar wind model opacities, and 
thereby illustrates the relative importance of two effects 
on the opacity: the elemental abundances, and the ion- 
ization balance. The solid line gives the same solar abun- 
dance wind model described in the previous paragraph, 
while the dashed and dotted lines give models particular 
to C Pup, using non-solar abundances specific to the star 
derived from detailed CMFGEN modeling (J.-C. Bouret 
et al., in preparation). The dashed line gives the ac- 
tual opacity at ~ 2R^ from this CMFGEN model, while 
the dotted fine uses the simplified ionization structure of 
the solar abundance wind opacity model. The fact that 
the realistic CMFGEN ( Pup model and the simplified 
version are so similar indicates that a simple ionization 
balance is typically adequate to describe wind opacity in 
many cases, as long as the dominant ionization stages 
are relatively accurate. On the other hand, the differ- 
ence between the ( Pup models and the solar abundance 
model shows that the opacity model depends strongly on 
the abundances of the most common elements other than 
H and He (typically C, N, and O). Note that the Bouret 
et al. abundances for ( Pup are both non-solar in the ra- 
tio of CNO and also have sub-solar metallicity; it is the 
sub-solar metallicity that accounts for the lower opacity 
of the C Pup models at short wavelengths compared to 
the solar abundance model. 

We have made one important simplification in our 
modeling: in Equation (j4]), and throughout this section, 
we have assumed that the opacity is independent of ra- 
dius. As shown in Figure [71 moderate changes in wind 
ionization do not strongly affect the opacity, so in most 
cases this is a justified approximation. The important ex- 
ception is the ionization of helium; in sufficiently dense 
winds, He++ may recombine to IIe+ in the outer part of 
the wind, which ^[[re atly increases the opacity, esp ecially 
at long wavelengths (|Pauldrachlll 98 71 : IHillier et al.l [l993). 
As long as the change in ionization occurs sufficiently far 
out in the wind, geometrical effects as described in Sec- 
tion [2] are not important, and the absorption due to He+ 
can be treated as an additional slab between the X-ray 
emitting regions and the observer, i.e., using T = e~^ . 

4. MODEL IMPLEMENTATION 
The numerical evaluation of Equation ([1^ 
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is not pro- 
hibitively expensive, but it is typically not fast enough 
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Figure 6. Comparison of a n eutral interstellar a bsorption model 
(including dust grains), tbabs (|Wilms et al.ll2000l '): and an O star 
wind model with an assumed simple ionization structure. The 
opacity ^c(A) is scaled with A^ to allow a better comparison of the 
two models at short w avelengths. Both models use solar abun- 
dances (Aspl und et al. ' 2009). Note that even below the 0^+ edge 
near 20 A, the wind opacity model is still about 40% lower than 
the ISM model, mainly due to the ionization of H and He in the 
wind model. 
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Figure 7. Comparison of three different model O star wind opac- 
ities: the solar abundance O star wind model (with assumed sim- 
ple ionization structure) shown in Figure [6] (solid); a CMFCEN 
model appropriate for (" Pup (dashed) using the abundances of J.- 
C. Bouret et al. (in preparation); and a simplified version of this 
CMFCEN model (dotted; red in the electronic version) where the 
ionization balance is the same as in the solar abundance model. 
Note that the Bouret et al. abundances for ^ Pup are subsolar as 
well as having altered CNO abundance ratios. 

to ahow its use in an automate d spectral fitt ing rou- 
tine, such as that in XS PEC12 (j ArnaudI 119961 ) or ISIS 
(|Houck fc Denicolall2QQQf ). It is thus preferable to com- 
pute the transmission on a grid in r* for a given set of 
wind parameters. 

Given a tabulation of the model wind opacity, as de- 
scribed in Section [3l in addition to the tabulation of 
T{r^)^ one may then calculate the transmission as a func- 
tion of wavelength, T(A), with only one free parameter, 
the characteristic wind mass column density S*. This 
parameter is analogous to the neutral hydrogen column 
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Figure 8. Transmission as a function of wavelength for ionized 
wind absorption model (windtabs, black) and for neutral slab ab- 
sorption (tbabs, gray; red in the electronic version). Three values of 
absorbing column are given; for windtabs, the degree of absorption 
is specified by the characteristic mass column density E*, while for 
tbabs it is given simply by the mass column E. 

density in a slab absorption model such as the XSPEC 
wabs or tbabs ISM absorption models (which are also 
sometimes used to approximate wind absorption). 

We have implemented this as a local model for XSPEC 
120. A user can calculate T{r^) for a given set of pa- 
rameters (i.e., /3, Rq), with the results stored in a FITS 
table. These implicit model parameters may be varied 
by computing additional tables of T{r^) for each set of 
parameter values. However, the absorption model is not 
very sensitive to these parameters over the range typi- 
cally inferred for winds of massive stars. The calcula- 
tion of T(r*) is controlled by a simple python script, and 
computation of a table for a given set of parameters can 
be accomplished in several seconds on a modern work- 
station. The model opacity must also be supplied as a 
FITS file; different model opacities may be swapped in 
at run time. When windtabs is used in a spectral fitting 
program, model transmission is calculated as a function 
of energy or wavelength using the supplied FITS tables 
and the one free model parameter, S*. 

The elemental abundances, which enter into the trans- 
mission through their effect on the opacity, cannot be 
varied as fit parameters in our model. This is a choice 
we have made in the model implementation, both for 
computational ease and simplicity of user interface, and 
because there is not enough information in X-ray spectra 
to constrain elemental abundances through modeling of 
absorption of X-ray emission alone. Abundances should 
be inferred by other means, e.g., from global fits to UV 
and optical spectra, or from fits to X-ray emission line 
strengths. These abundances can be used to compute a 
new opacity table for a given star. 

Figure [8] gives the model transmission for windtabs us- 
ing three different values of S*, and using our standard 
O star wind opacity model. For comparison, this figure 
also shows the transmission for the neutral absorption 

^ Source code is freely available under the 
General Public License, and it may be ob- 



tained 



alon^ 



with 



installation 



instructions 



at 



|http : //heasarc . gsf c . nasa . gov/docs/xanadu/xspec/models/windprof . 



model tbabs for comparable mass column densities S. 
Note that S refers simply to a slab mass column density, 
while S* refers to a characteristic mass column density 
in the context of a stellar wind (see Equation (fTTj) ). 

5. DISCUSSION 

5.1. Advantages of windtabs over exospheric or 
exponential attenuation 

The exact method we have presented here for modeling 
the X-ray radiative transfer for a distributed source em- 
bedded within a stellar wind is crucial for analyzing the 
X-ray emission observed from O stars. Compared to the 
exponential, neutral slab absorption (excess over ISM) 
model that is usually employed, the windtabs model ac- 
curately reproduces the much more gradual decrease in 
transmission with increasing wind column density and 
opacity. This can be seen in Figure [H where the exponen- 
tial transmission model shows an unrealistically sudden 
decrease in transmission when the wind becomes opti- 
cally thick. 

Because the opacity of the bulk wind is a relatively 
strong function of wavelength, the inaccuracy of the ex- 
ponential transmission model will lead to errors in the 
broadband spectral energy distribution of a model ap- 
plied to individual O stars, leading to misinterpretations 
of the associated emission model components. This ap- 
pears to be the case in the stud y of Chandra grating 
spectra in ' Zhekov fc Pa lla' (2007) , which invokes excess 
exponential, neutral ISM absorption to account for the 
assumed wind attenuation. Most likely because the ex- 
ponential transmission radiation transport model, and 
also the neutral ISM opacity model, significantly overes- 
timates the degree of attenuation, these authors find ab- 
sorption beyond the ISM column for only one star, even 
though some additional wind attenuation is expected for 
most of the stars in their sample. Additionally, although 
the authors do not comment on it, their determination of 
elemental abundances for each of the stars shows a consis- 
tent trend of abundance correlated with the wavelength 
of the dominant line or lines from each element. Such 
an effect would be expected if the wavelength-dependent 
wind attenuation is not accurately accounted for. 

It has long been noted that the exponential attenua- 
tion treat ment is no t well suited to modeling OB star 
X-rays. iHillier et al.l (|1993f ) used an exact treatment in 
their modeling of ( Pup, discarding the possibility of a 
coronal model on the basis of the strong soft observed 
X-ray flux of C Pup, toget h er wi th its relatively high 
mass- loss rate. I Cohen et al.l (|1996f ) showed that an exo- 
spheric treatment, rather than an exponential treatment, 
is important for understanding the observed EUV and 
soft X-ray emission from the early B giant, e CMa. The 
exosp heric approximation was used by Qwocki fc CohenI 
(|1999f ) to model the effect of wind attenuation of X-rays 
in order to explain the observed La^/Lboi ^ 10~^ rela- 
tionship and its breakdown in the early B spectral range 
where hot star winds become optically thin to X-rays. 
An ex ospheric treatment was also used bv lOskinova et aP 
(|2001[ ) to analyze the variability of X-rays from optically 
thick WR winds. And the exospheric framework forms 
the basis for the "optical depth unity" relationship for X- 
rays in O stars, where the forbidden-to-intercombination 
Mig^ ratios of helium-like ions are claimed to imply for- 
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mation radii that track the optic al depth unity radius 
as a f unction of X-ray wavelength (IWaldron fc Cassinellil 
[2003). 

However, the exospheric treatment underestimates the 
attenuation of the wind. If data need to be analyzed with 
a high degree of accuracy, then a more realistic treatment 
of X-ray radiation transport through the wind must be 
used; one that takes the inherently non-spherically sym- 
metric nature of the problem into account. This is espe- 
cially true when the location of the X-ray emission emerg- 
ing from the wind is important, as with the interpretation 
oijji ratios. As Figure [2] shows, the relative contribution 
of the emergent X-ray flux from different wind regions is 
incorrectly estimated in the exospheric approach, and as 
Figure m shows, the emergent flux has significant contri- 
butions from a wide range of radii. In fact, a number 
of previous works have implemen ted accurate X-ray ra- 
diative transfer prescriptions (e.^.JPaulc rach et aIll2QQll : 
lOwocki fc CohenI l2QQll : 10 skinova et al .1120061 ): however, 
none of these works provide an implementation that is 
available for use in other contexts. 

The windtabs model we have introduced here is not 
only more accurate in terms of the radiation transport 
than slab or exopheric treatments, but it has two addi- 
tional advantages that recommend its adoption for rou- 
tine X-ray data analysis and modeling of O stars. First, 
it is easy to use and has only a single free parameter, 
the characteristic mass column density S*, from which a 
mass-loss rate can be readily extracted. And second, it 
incorporates a default wind opacity model that is signifi- 
cantly different from, and much more accurate than, the 
neutral ISM opacity models that are usually used. Ad- 
ditionally, alternate user-calculated opacity models are 
easy to incorporate. 

5.2. Implications for trends in X-ray hardness with 
spectral type 

One application of windtabs to the interpretation of X- 
ray spectral data is for the analysis of the X-ray spectral 
hardness trend versus optical spectral s ubtype recently 
noted in Chandra grating spectra by Walborn et al. 
( l2009l ). They report "the progressive weakening of the 
higher ionization relative to the lower ionization X-ray 
lines with advancing spectral type, and the similarly de- 
creasing intensity ratios of the H-like to He-like lines of 
the a ions." The correlation of the overall X-ray spec- 
tral hardness with spectral subtype described by Wal- 
born et al. appear to be a direct result of wavelength 
dependent absorption effects. A detailed analysis is in 
preparation, but here we show a single suite of models 
in which the only variable parameter is the characteris- 
tic wind column density, S*. A single emission model, 
combined with windtabs attenuation, reproduces the ob- 
served broadband trend quite well. 

Figure [9] shows X-ray grating spectra for seven O giants 
and supergiants in the left-hand column, with the earliest 
spectral subtype (03. 5) on the top, and the latest (BO) 
on the bottom, as in IWalborn et al] (j2009f ). The later 
spectral subtypes clearly have more soft X-ray emission, 
although the earlier subtypes still have non-negligible 
long- wavelength (A > 15 A) emission. In the middle 
column we show a four-temperature apec (jSmith et al.l 
l2001f ) thermal equilibrium emission model. We have cho- 
sen kT = 0.1, 0.2, 0.4, and 0.6 keV; the first three compo- 





Table 1 










Adopted Stellar Parameters 






Star 


Spectral Type^ 


A^H^ 


S*^ 


A^*^ 


HD 150136 


03.5 If* + 06V 


0.76 


0.073 


3.6 


CPup 


04 I(n)f 


0.01 


0.160 


7.2 


CPer 


07.5 III(n)((f)) 


0.115 


0.017 


0.76 


r CMa 


09 II 


0.056 


0.013 


0.58 


6 Ori 


09.5 II + B0.5 III 


0.015 


0.011 


0.49 


C Ori 


09.7 lb 


0.03 


0.019 


0.85 


e Ori 


BO la 


0.03 


0.020 


0.89 



^ Spectral types adopted by Wal born et al.l (|2009I V 
^ Interstellar medium column density (10"^"^ cm"-^). 
^ Characteristic wind mass column density (g cm~^). 
^ Characteristic wind number column density equivalent 
toE* (1022 cm-2). 

nents having equal emission measures and the hottest one 
having half the emission measure of the others. The same 
APEC model is used for all seven stars (with variable over- 
all normalization) and is multiplied by a windtabs model 
and a tbabs model (for neutral ISM attenuation). The 
column density of the t babs model is fixed a t the inter- 
stellar value taken from iFruscione et al.l ("1994), with the 
exception of HD 150136, for which we inferred the ISM 
column density from E(B-V) ()Maiz-ADellaniz et all2QQ4l : 
iMartins fc Ple3 I2QQ6I : IVuong et al.l l2QQ3f ) . Thediarac- 
teristic mass column density S* in windtabs is fixed at 
a value computed from the "cooking formul a'' the oret- 
ical mass- loss rate computed bv IVink et al.j ( 2001 ). us- 
ing the measured terminal velo city of iHaserl ()1995D and 
the modeled radii of .Martins et al.l (|2005f ). The stan- 
dard solar abundance wind opacity model (solid line in 
Figure [7j) was used in windtabs, and the APEC model 
abundances were set to solar. There are no free param- 
eters in these models, and the temperature distribution 
has not even been significantly optimized to match the 
data. The adopted parameters are listed in Table [TJ 

As the middle column of Figure [9] shows, the sim- 
ple, universal emission model reproduces the broadband 
trend very well. Trends in individual line ratios generally 
cannot be reproduced only by ac counting for the vary - 
ing attenuation, as pointed out bv IWalborn et all (|2009f ). 
but note that the Ne IX (13.5 A) to Ne X (12.1 A) ra- 
tio does indeed vary due only to differential attenuation 
among the earliest spectral subtypes. The right-hand 
column in Figure [9] shows the same emission and ISM 
attenuation models as in the middle column, but with 
excess exponential (neutral ISM) attenuation accounting 
for the wind absorption, again according to the wind col- 
umn densities expected from the adopted mass- loss rates, 
radii, and terminal velocities. The exponential attenua- 
tion trend seen in the right-hand column is too strong for 
the earliest spectral subtypes and too weak for the latest 
ones, where the different ISM column densities actually 
dominate the trend. 

The contrast between the windtabs and exponential 
models is quite stark, and indicates that the more re- 
alistic models should generally be used when analyzing 
X-ray spectra, both high-resolution and broadband. It is 
also impressive how much of the observed spectral hard- 
ness trend is explained by wind attenuation, in the con- 
text of a realistic model. Not only do quantitative analy- 
ses of the suggested line ratio trends have to be evaluated, 
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Figure 9. Left column: a sequence of Chandra spectra of O giants and supergiants from IWalborn et al.l (|2009l ): middle column: multi- 
temperature thermal emission model with windtahs wind absorption model; right column: same model as middle column, but with thahs 
neutral slab absorption model. 



but a global spectral modeling that allows for both emis- 
sion temperature variations and wind attenuation vari- 
ations should be undertaken in order to disentangle the 
relative contributions of trends in emission and absorp- 
tion to the overall, observed trend in the spectral energy 
distributions. 

6. CONCLUSIONS 

We have presented an exact solution to the radiation 
transport of X-rays through a spherically symmetric, 
partially optically thick O star wind, and shown that 
it differs significantly from the commonly used slab ab- 
sorption and exospheric models. Specifically, the trans- 
mission falls off much more gradually as a function of 
fiducial optical depth in the windtahs model as compared 
to the exponential model, leading to more accurate as- 
sessments of wind column densities and mass-loss rates 
from fitting X-ray spectra. As one example of the utility 
of windtahs^ we have shown that when this more accu- 
rate model is employed, differential wind absorption can 
explain most of the observed trend in OB star X-ray 
spectral hardness with spectral subtype, and even may 
explain some of the line ratio trend. 

The windtahs model has been implemented as a custom 
model in XSPEC, and is as easy to use as the various ISM 
absorption models, having only one free parameter. In 
addition to the significantly improved accuracy of the ra- 



diation transport, windtahs has several other advantages. 
It incorporates a default opacity model much more ap- 
propriate to stellar winds than the neutral element opac- 
ity model used in ISM attenuation codes. Users can eas- 
ily substitute their own custom-computed opacity mod- 
els. And the fitted mass column density parameter for 
windtahs allows for the user to extract a mass-loss rate 
directly from their fitting of X-ray spectra of OB stars. 
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